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ABSTRACT 



An algorithm and its first implementation in C# are presented for assembling arbitrary quantum circuits on the 
base of Hadamard and Toffoli gates and for constructing multivariate polynomial systems over the finite field 
Z 2 arising when applying the Feynman's sum-over-paths approach to quantum circuits. The matrix elements 
determined by a circuit can be computed by counting the number of common roots in Z-i for the polynomial 
system associated with the circuit. To determine the number of solutions in Z 2 for the output polynomial system, 
one can use the Grobner bases method and the relevant algorithms for computing Grobner bases. 
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When applying the famous Feynman's sum-over-paths approach to a quantum circuit, as it was shown in, 1 one 
can associate with an arbitrary quantum circuit a polynomial system over the finite field Z 2 . For N input qubits 
the system contains N + 1 polynomials: each of the input qubits adds a certain polynomial to the system, 
and one more polynomial arises from the phase of a classical path through the circuit. Thereby, elements of 
the unitary matrix determined by the quantum circuit under consideration can be computed by counting the 
number of common roots in Z2 for the polynomial system associated with the circuit. Given a quantum circuit, 
the polynomial set is uniquely constructed. 

In this paper which is a substantially extended version of paper 2 we describe an algorithm for building an 
arbitrary quantum circuit made of Hadamard and Toffoli gates and for constructing the associated polyno- 
mial system. The algorithm implemented in the form of C# 6 program called QuPol (abbreviation of Quantum 
Poly nomials). The program provides a user- friendly graphical interface for building quantum circuits and gen- 
erating the associated polynomial systems. 

To build an arbitrary quantum circuit, it is represented as a rectangular table with N x M cells, where N 
the number of the input qubits, M the number of the circuit cascades plus one - for the input iV-qubit vector. 
Each cell contains an elementary gate taken from the elementary gate set. The last set of gates consists of wires 
(identities) , arithmetic operations (multiplication and addition) and the Hadamard operation. A user of QuPol 
can construct a quantum circuit by selecting required elementary gates from the menu bar of the program and by 
placing them in the appropriate cells. As a result, the constructed table provides a visual representation of the 
circuit. A recursive evaluation procedure built-in QuPol operates on the table column- by-column and generates 
N polynomials (one for each row of the table) as well as the phase polynomial in n variables, where n the number 
of Hadamard gates in the circuit. 

To count the number of solutions in Z 2 for the output polynomial system, one can convert it into an appro- 
priate Grobner basis form which has the same number of solutions. Being invented 40 years ago, 3 the Grobner 
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1. INTRODUCTION 



bases has become the most universal algorithmic method for analyzing and solving systems of polynomial equa- 
tions. 4 In particular, construction of the lexicographical Grobner basis substantially alleviates the problem of 
the root finding for polynomial systems. To construct Grobner bases one can use, for example, efficient involutive 
algorithms developed in. 5 Our QuPol program together with a Grobner basis software provides a tool to analyze 
quantum computation by applying methods of modern computational commutative algebra. We illustrate this 
tool by example from. 1 

The structure of the paper is as follows. In Section 2 we outline shortly the circuit model of quantum 
computation. Section 3 presents the famous Feynman's sum-over-paths method applied to quantum circuits. 
In Section 4 we consider a circuit decomposition in terms of the elementary gates and show how to build an 
arbitrary circuit composed from the Hadamard and Toffoli gates. These two gates form a universal gate basis. 7 ' 8 
Some features of the algorithm implementation are briefly described. Section 5 demonstrates a simple example of 
handling the polynomials associated with a quantum circuit by constructing their Grobner basis and computing 
the circuit matrix. We conclude in Section 6. 

2. QUANTUM CIRCUITS 

To compute a reversible Boolean vector-function / : 1% — > Z£ , one applies an appropriate unitary transformation 
Uf to an input state |a) composed of of n qubits 9 : 

|b) - U f |a) , |a) , |b) e C 2 " . 

The output state |b) is generally not required result of the computation until somebody observes it. After that 
the output state becomes classical and can be used anywhere. 

Some elementary unitary transformations are called quantum gates. A quantum gate acts only on a few 
qubits, on the remaining qubits it acts as the identity. One can build a quantum circuit by appropriately aligning 
quantum gates. In so doing, the unitary transformation of the circuit is the composition of its elementary unitary 
transformations: 

U f = U m U m -i---U 2 U 1 . (1) 

A quantum gate basis is a set of universal quantum gates, i.e. any unitary transformation can be presented 
as a composition of the gates of the basis. As well as in the classical case, there can be different choice of the 
basis quantum gates. 9 For our work it is convenient to choose the particular universal gate basis consisting of 
Hadamard and Toffoli gates. 7, 8 

The Hadamard gate is a one-qubit gate. It turns the computational basis states into the equally weighted 
superpositions 

:|0)~i=(|0> + |l», 
# :|1>~^(|0>-|1» 

which differ in the phase factor at |1). 

The Toffoli gate is a tree-qubit gate. The input bits x and y control the operation on bit z, and the Toffoli 
gate acts on the computational basis states as 1 

(x,y,z) i ^ (x,y,z@xy) . 

An action of a quantum circuit can be described by a square unitary matrix whose entry (b| Uf |a) yields the 
probability amplitude for the transition from the initial quantum state |a) to the final quantum state |b). The 
matrix element is represented in accordance to the gate decomposition (1) of the circuit unitary transformation 
and can be calculated as sum over all the intermediate states a i; i = 1,2, . . . m - 1: 



(b| U f | a) - ]T (b| U m |a m _!) ■ • • < ai | Ux |a> 



3. SUM-OVER-PATHS FOR QUANTUM CIRCUITS 



To apply the famous quantum-mechanical Feynman's sum-over-paths approach to calculation of the matrix 
elements for a quantum circuit, 1 we replace each quantum gate in the circuit under consideration by its classical 
counterpart. The trick here is to determine the corresponding classical gate for the quantum Hadamard gate 
since its action at any input value or 1 must generate or 1 with the equal probability. To take this into 
account, the output of the classical Hadamard gate can be characterized by the path variable i £ Z 2 .' Its value 
determines one of the two possible paths of computation. Thereby, the classical Hadamard gate as 



Oi 



aj, x £ Z 2 , 



and the classical Toffoli gate acts as 



(oi, o 2 , a 3 ) i ► (ax, a 2 , a 3 © aia 2 ) 



where © denotes addition modulo 2. 

Fig. 1 shows an example of quantum circuit (taken from 1 ) and its classical counterpart. The path variables 
Xi comprise the (vector) path x = (x\, x 2 , Xs, X4) T <G Zf. 
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Figure 1. From quantum to classical circuit 



A classical path is defined by a sequence of classical bit strings a, ai, a 2 , . . . , a m = b obtained from application 
of the classical gates. Each set of values of the path variables Xi gives a sequence of classical bit strings which is 
called an admissible classical path (Fig. 2). Each admissible classical path has a phase which is determined by 
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Figure 2. Admissible classical paths 

the Hadamard gates applied. 1 The phase is changed only when the input and output of the Hadamard gate are 
simultaneously equal to 1, and this derives the formula 

(f(x) = ^ input • output . 

Hadamard gates 

with the sum evaluated in Z 2 . As to Toffoli gates, they do not change the phase. 



In our example the phase of the path x is 

ip(x) = a 1 x 1 a 2 x 2 ® Xxx 3 x 4 (a 3 x x x 2 ) ■ 



The matrix element of a quantum circuit is the sum over all the allowed paths from the classical state a to b 

x:b(x)b 

where h is the number of Hadamard gates. The terms in the sum have the same absolute value but may vary in 
sign. 

Let No be the number of positive terms in the sum and N\ be the number of negative terms: 



N = | { x | b(x) = b A p(x) - } 
JVi = j { x | b(x) = b A ip{x) = 1 } 



(2) 
(3) 



Hence, Nq and Ni count the number of solutions for the indicated systems of n + 1 polynomials in h variables 
over Z2. Then the matrix element may be written as the difference 



(b| U f |a> 



1 



(N - JVi) . 



(4) 



4. ALGORITHM AND PROGRAM 

For building arbitrary quantum circuits from Hadamard and Toffoli gates, we use the set of elementary gates 



+ A v A V A V 

E = {I,I,I,I,M,M,A,A,H} 
shown on Fig. 3. Furthermore, we represent a circuit as a rectangular table (Fig. 4, left). 
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Figure 3. Elementary gates 

Each cell in the table contains an elementary gate so that the output for each row is determined by the com- 
position of the elementary gates in the row (taking into account recursive dependences). Hence, each elementary 
unitary transformation Uj is represented by an n-tuple of elementary gates. 

Fig. 3 shows action of the elementary gates from the set (5): the identities, the multiplications, the additions 
modulo 2, and the classical Hadamard gate. The identity I just reproduces its input. The identity-cross 
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We want to assemble a circuit 
with 3 rows and 4 columns 
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Figure 4. Decomposition into elementary gates (left) and building a sample circuit, step (right) 
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Figure 5. Building circuit, steps 1 (left) and 2 (right) 



+ 

I reproduces also its vertical input from the top elementary gate to the bottom one and vice-versa. Every 

V A A 

identity-down / and identity-up / have two outputs - horizontal and vertical. The multiplication-up M and 

v 

multiplication-down M perform multiplication of their horizontal and the corresponding vertical inputs. The 

A V 

addition-up A and addition-down A are applied similarly. Each Hadamard gate outputs an independent path 
variable irrespective of its input and can give a nonzero contribution to the phase. 

To build a circuit, we define, first, an empty table of the required size (Fig. 4, right). At this initialization 
step the output qubit values and phase are not determined. As the second step we place the required elementary 
gates in appropriate cells of the first column and construct the circuit polynomials (Fig. 5, left). Now the output 
is obtained by applying the inserted elementary gates to the input. The phase is also calculated. Then we 
proceed the same way with the second column (Fig. 5, right), then with the third and forth columnns (Fig. 6). 

A circuit is represented in our program QuPol as two 2d- arrays: one for the elementary gates and another 
for the related polynomials. The phase polynomial is separately represented. The following piece of code 
demonstrates construction of the circuit polynomials 



for each Column in Table of Gates 
for each Gate in Column { 
construct Gate Polynomial; 
if Gate id Hadamard 

reconstruct Phase Polynomial; } 
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Figure 6. Building circuit, steps 3 (left) and 4 (right) 

It should be noticed that Figs. 5 and 6 show a shortened version of the elementary gates composition. Actually, 
the final polynomials look like: 

h = A(HlHa 1 ,M(lM(Ha 2 ,lHa 1 ),lHA{Ia 3 ,M(Ha 2 ,lHa 1 )))) =x 3 ® x 2 x 4 , 
b 2 = M{I M(Ha 2 ,} H ai ),I HA(Ia 3 ,M(Ha 2 J H ai ))) = x 2 , 
h = lHA(Ia 3 ,M(Ha 2 jH ai )) = x 4 . 

The method for constructing a gate polynomial is essentially recursive because of the need to go up or down 
for some gates. 

Any circuit is saved in the program as two files. One file is binary and contains the circuit itself. Another 
file has a text format. It contains the circuit polynomials in a symbolic form. The program allows one to save 
polynomials in several formats suitable as inputs for computer algebra systems such as Maple or Mathematica 
for the further processing. It is also possible to load the saved circuit back into the computer memory for its 
modification. 

The part of our code for the name space PolynomialJVIodulo_2 written in C# can be also used indepen- 



Class Polynomial: Enumerable 



• void AddTo (Polynomial polynomial) 

• void MultiplyBy (Polynomial polynomial) 

• override string ToString () 

• bool FromString (string str) 

• float AsGraphics (...) 

• void Collect () - collect terms 

• int Count () - get number of monomials 



Figure 7. Class Polynomial 



Methods 




dcntly on our program. This name space contains a set of classes for handling polynomials over the field Z2. 
Class Polynomial (Fig. 7) is a list of monomials, class Monomial (Fig. 8) is a list of letters, class Letter (Fig. 9) 



Class Monomial : Enumerable 



■ void AddTo (Monomial monomial) 

• void MultiplyBy (Monomial monomial) 

• void MultiplyBy (Letter letter) 

• override string ToString () 

• bool FromString (string str) 

■ float AsGraphics (...) 

■ int Count () - get number of letters 



Methods 



Figure 8. Class Monomial 
is an indexed letter provided with a positive integer superscript (power degree). 



Class Letter: IComparable 



• char letter 


Properties 


• int index 




• int power 





■ override bool Equals (object obj) 

■ int CompareTo (object obj) 
override string ToString () 

■ bool FromString (string str) 

■ float AsGraphics (...) 



Methods 



Figure 9. Class Letter 



5. CALCULATING CIRCUIT MATRIX ELEMENTS 

A system generated by the program is a finite set F C R of polynomials in the ring 

R = Z 2 [a i ,b j ][x 1 , ...,x h ] , a i ,b j eZ 2 , i,j = l,...,n 



(6) 



with h polynomial variables and 2n binary coefficients. Let iVo and N\ denote, as in (2-3), the number of roots 
in Z 2 of the polynomial sets F and F\ , respectively, where 

F = {/,..., f k , <p}, F 1 ={/,... ,/fc^el}. (7) 



Then the circuit matrix is given as 

(b\U\ a ) = -±= n {N -N 1 ), 

where h is the number of Hadamard gates in the circuit. 

To compute N n and Ni one can convert F and F\ into an appropriate triangular form 10 providing elimination 
of the path variables x\, ...,Xh- One of such triangular forms is the pure lexicographical Grobner basis that can 
be computed, for instance, by means of the Buchberger algorithm, 3,4 by the Faugere algirthms 11 ' 12 or by our 
involutive algorithm. 5 

For the circuit shown on Fig. 1 we obtain the following polynomial set: 

fi = x 2 x 4 © x 3 © bi , 
h = x 2 © b 2 , 
h = x 4 © b 3 , 

ip = x\x 2 © xix 3 © aiXi © a 2 x 2 © 0,3X4 



(8) 



The lexicographical Grobner bases for polynomial systems F and F\ and for the ordering x\ >~ x 2 >~ x 3 >- x 4 
on the variables are given by 



G : ^ 



9i = (01 © &i)a;i © a 2 b 2 © a 3 & 3 , 

.92 = x 2 © 6 2 , 
53 = x 3 © 61 © 6263 , 
k .93 = X 4 © & 3 , 



Gi: <^ 



9i = (ai © © a 2 6 2 © 03^3 © 1 

.92 = 2:2 © b 2 , 

33 = ^3 ffi &1 ffi &2&3 , 

, 33 = x 4 © h . 



(9) 



These lexicographical Grobner bases immediately yield the following conditions on the parameters: 



G : a\ © 61 = a 2 & 2 © 03^3 = , 
Gi : ai © 61 = A a 2 b 2 © a 3 b 3 



1. 



(10) 
(11) 



It is immediately follows that if conditions (10) are satisfied then the polynomial system Go (resp. Fq) has 
two common roots in Z 2 and G\ (resp. F\) has no common roots, and, vise-versa, if conditions (10) are satisfied 
then Go has no roots and Gi has two roots. In all other cases there is one root of Go and one root of G\. 

In that way, the 8x8 matrix for the circuit of Fig.l is easily determined by formulae (4) with 7V and N\ 
defined from systems (9). Table 1 explicitly shows entries of this matrix when its rows and columns are indexed 
by the values of the input bits 01, a 2 , ft3 and output bits b\,b 2 ,b 3 , respectively. 



6. CONCLUSION 

We presented an algorithm for building an arbitrary quantum circuit composed of Hadamard and Toffoli gates 
and for constructing the corresponding polynomial equation systems over Z 2 . The number of common roots Z 2 
of polynomials in the system uniquely determines the circuit matrix. The algorithm has been implemented in as 
a C#-program QuPol. 

To count the number of polynomial roots one can use the universal algorithmic Grobner basis approach. In 
the framework of this approach the polynomial system under consideration is converted into a triangular form 
that is convenient for computing the number of roots. 
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Table 1. The circuit matrix for polynomial system (8) 

General theoretical complexity bound for computing Grobner bases (see, for instance, 13 ) is double exponential 
in the number of polynomial variables Xi (the number of path variables, i.e. Hadamard gates, in our case). 
However, the due to the fact that all the path variables Xi, the input and output parameters aj,bk as well as 
the numerical coefficients of the polynomials are elements in the finite field Z 2 the Grobner basis computation 
is sharply simplified. Thus, based on the algorithm in, 12 the analysis of costs for computing Grobner bases for 
polynomials over Z 2 , arising in cryptography, revealed 14 only a single exponential complexity in the number of 
equations. 

Therefore, the above presented algorithm together with an appropriate software for computing a Grobner 
basis and for counting its number of common roots in Z 2 provides a tool for simulating quantum logical circuits 
and for estimation of their computational power. For this purpose we are planning to create a special module 
in the open source software GINV. 16 At present, GINV contains implementation in C++ of our involutive 
algorithm 5 for computing Grobner bases over the field of rational numbers and over the field Z p where p is 
a prime number. By this reason, the current version of GINV can be used for computing Grobner bases for 
polynomial systems (7) only if one considers 2n parameters a,i,bj in (6) as variables. In this case, to compute 
a triangular Grobner basis form of the circuit polynomial, one should use an elimination term order for those 
parameters such that any term containing path variables Xk is higher than that any term containing parameters 
only. 

Such computational scheme, however, increases the number of polynomial (path) variables by 2n what may 
lead to exponential slowing down, as was pointed to above. That is why we expect that an adaptation of the 
built-in algorithms and data structures to the polynomial ring (6) with parametric coefficients will substantially 
increase computational efficiency of GINV in its application to analysis and simulation of quantum circuits. 

Recently another algorithmic approach was suggested 15 to work with polynomial systems over finite fields. 
It seems interesting to investigate computational efficiency of the approach in 15 in comparison with that based 
on Grobner bases. 
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